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The problem of intermittency in developed hydrodynamic 
turbulence is considered. Explicit formulae taking into ac- 
count effects of finite size of the inertial range are presented 
for the whole set of intermittency exponents. The formulae 
fit pretty well experimental data whose apparent discrepan- 
cies are attributed to different sizes of the inertial ranges in 
different experiments. Further predictions are given that can 
be verified by already existing experimental technique. 

• INTRODUCTION: Developed hydrodynamic 
turbulence was studied by many scholars of our century. 
The subject is complicated and a theory that fits a broad 
experimental database, having clear physical background 
and real predictive power is still wanted. The problem is 
known to be closely connected with the problem of blow- 
up of velocity gradients in the Euler fluid. The latter 
problem is also complicated and unsolved yet. Neverthe- 
less, the joint study of these problems can be helpful and 
based on it one can arrive to a new qualitative picture of 
the intermittency and also to new quantitative results. 

• INTERMITTENCY EXPONENTS: Let e{r) 
be the local density of energy dissipation rate in a 
strongly turbulent fluid. The time argument t is omitted 
here and below as all quantities considered at the same 
time t and as the forthcoming space averaging removes 
all time dependence. The energy dissipation rate can be 
smoothed over a spatial scale x by means of a "filter func- 
tion" F(p): e(r,x) = f dpF(p)e(r + xp), f dpF(p) = 1. 
The result depends on the filter F. This dependence is 
important and will be described elsewhere. For now, IE- 
dependence is considered assuming that the filter satisfies 
some basic requirements. For x inside the inertial range, 
i.e., for L x h, where L is the energy-containing 
and h is dissipative (Kolmogorov) scale of the turbulence, 
some scale-invariance properties are expected. In partic- 
ular, it is assumed that for x and y both inside the inertial 
range 



(( £ (r,y)/e(r,x)) q ) = (x/y) 



(1) 



where the angle brackets signify space averaging. The 
scale invariance hypothesis is that the exponents p q (x/y) 
depend just on the ratio x/y rather than on both vari- 
ables x and y. These exponents are usually referred 
to as the intermittency exponents, because they all are 
zero in the Kolmogorov model of turbulence [1] that as- 
sumes spatially uniform energy dissipation e. There are 
several definitions of the intermittency exponents in the 
literature, which sometimes assumed to be equivalent, 



but, strictly speaking, should be distinguished from each 
other. 

It is convenient for further applications to change vari- 
able so that x/y = 2" and p q (x/y) = p n ,q- As n — ► oo 
(which limit makes sense just for infinitely broad iner- 
tial range), the intermittency exponents become scale- 
independent: q — ► n q . There are many experiment- 
based evaluations of p q with widely scattered results (and 
there are many phenomenological models trying to fit 
these evaluations). This is basically due to effects of fi- 
nite size of the inertial range, which are important in all 
realistic situations. The importance of finite size correc- 
tions to the intermittency exponents strongly complicates 
the problem. It is possible however to introduce physical 
quantities which are less sensitive to the finite size effects 
than intermittency exponents. 

• FRACTAL DIMENSIONS: Consider the spa- 
tial domain H(a, x, y) defined by condition 

e(f,y)/e(f,x) > (x/yf . (2) 

Its average volume per cube x 3 depends on y < x as 

V(a,x,y) = y 3 (x/yf^M = x 3 (x/yf^M- 3 , (3) 

which employs the scale invariance hypothesis again. The 
quantity V/y 3 can be interpreted as the smallest average 
number of cubes y 3 per cube x 3 required to cover domain 
H(a,x,y). Then, D(a,x/y) = D n (a) is the dimension- 
ality of H(a, x, y). It cannot exceed 3 and it is a decreas- 
ing function of a. There is a„ o, such that D n (a) = 3 
for a < a n fi and D n (a) < 3 for a > a„ o- There is also 
®n,M, such that D n (a) > for a < a n ,M and D n (a) < 
for a > a n ,M- It follows from (l)-(3) that, for a positive 
<1, 



Hn, q > aq + D n (a) - 3 



(4) 



The difference between left-hand side of (4) and its right- 
hand side is a non-negative function of q (and other pa- 
rameters). It has a minimum over q. Provided a is such 
that the minimum is reached at a finite positive q, its 
position is determined by p! n = a (the prime signifies 
(/-derivative) and its value is zero, so that there 



D n (a) = 3 + fi„ 



(5) 



At n — ► oo the dimensionality D n q 
independent limit D q . There are other 



tends to its scale- 
definitions of 



fractal dimensions D q in the literature which should be 
distinguished from that given above. 
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Formula (5) expresses D nq in the terms of n n>q . When 
D nq is known, jj, n q can be restored by 



l^"n,q 



i(3 - D nAl )jq\ + ct nfi q 



(6) 



that takes into account condition = 0. The integral 
in (6) converges at qi — ► 0, as D„ o = 3 (see (5) at q=0) 
and it is the maximum value of D n q . It also makes clear 
that the above a n o corresponds to q = 0. Another useful 
form of (6) is 



l^"n,q 



L (3 - D nAl )/ql + Hn,i Q 



(7) 



The first intermittency exponent is usually assumed to 
be zero, but in the framework of given definitions there 
is no reason for conclusion i = at all n, although it 
is correct at n = oo. 

• HOT SPOTS: First consider the limit of infinitely 
broad inertial range. As n — » oo the quantities D n q, l^n,q 
and a n q = pt n tend to their extreme scale-independent 
values D q , pt q and a q respectively. At fixed q, the quan- 
tity N = V/y 3 (defined by (3) and interpreted after it) 
depends of n like N = 2 nDq . 

There is a physical interpretation of such a depen- 
dence. Recall that N is the number of sub-cubes with 
side y = x/2" per cube x 3 , required to cover domain 
H(a,x,y) where (2) is satisfied. Condition (2) selects 
the sub-cubes in which the energy dissipation rate ex- 
ceeds certain level, i.e. the "hot spots" of certain in- 
tensity. One can assume that, as n increases, the hot 
spots split into those of smaller size which then do the 
same at larger n producing finally a fractal structure. 
This picture agrees qualitatively with the clustering of 
energy dissipation in turbulence. There is also a natu- 
ral dynamic process underlying such a picture. It is a 
branching blow-up process — fractal collapse of velocity 
gradients in Euler fluid. If branching happens on aver- 
age after each increase n by some A, then it should be 
N = 2 n / A branches (i.e. hot spots) at given n. The com- 
parison with the previous expression for the number N 
shows, that D q = 1/A can be interpreted as the branch- 
ing coefficient for fractal collapse. 

Convexity of n n>q as the function of q implies that 
a n q = /i' n is an increasing function of q. Since D n (a) 
is a decreasing function of a and D n (a n q ) = D n q , one 
concludes that D n q is a decreasing function of q. The 
physical interpretation is that the branching delays col- 
lapse and reduces the strength of singularities. Taking 
higher moments q selects more singular, i.e. less branched 
structures. At sufficiently large q (and n — ► oo), just 
perfect non-branched structures are selected that corre- 
spond to the fastest possible blow-up of velocity gradi- 
ents in Euler fluid. For larger q (q > qm), one has 
D q = , n q = um 1 + 3 . The largest singularity expo- 
nent um can be determined from a nonlinear eigenvalue 



problem for self-similar blow-up solutions of Euler equa- 
tions. This is not done yet because of severe technical 
difficulties. Moreover, even the dimensionality of singu- 
larities is under discussion. The above picture assumes 
that singularities are point ones. There are also phe- 
nomenological models of turbulence (cited below) which 
assume line or plane singularities. 

• QUANTITATIVE PREDICTIONS: Provided 
certain universality properties of collapse branching, se- 
lection of more singular structures by taking higher mo- 
ments q should be a uniform process in q. Specifically, 
increasing q by some value should cause a (/-independent 
decrease in D q , i.e., it should be D' q = —7 = const, if the 
dimensionality is determined just by collapse branching. 
There is however a background that goes down to com- 
pensate for the increase of e in hot spots and secure the 
energy flux conservation. Indeed, it is easy to see that 
a = n' < 0. The background contribution to the di- 
mensionality is important at small q. The dimensional- 
ity (/-dependence at small q is quadratic in q, since q = 
is the point of D q maximum (Do = 3). This quadratic 
decrease in q is turned into the linear one at higher q 
where behavior of D q is already determined by the col- 
lapse branching. A simple formula for D q having such 
properties is 



D„ 



(8) 



This formula is valid for q < qm, while D q = for q > 
9m, 



Qm = 3/7 



1/2+ (l/4 + (a 7 /3) 2 ) 



i/2- 



1/2 



Substitution (8) in (7) at n = 00 gives 

Hq = 7 q [arsh (q/a) — arsh (1/a)] 



(9) 



(10) 



where arsh^ = ln(£ + \/ £, 2 + 1) is the hyperbolic arcsin. 
There is a connection between parameters a and 7 (which 
is specified below), so that all quantities at n = 00 can 
be expressed in the terms of a single parameter, say 7. 

Having \i q known, one can proceed with the calcula- 
tion of finite size corrections. Their calculation, taking 
into account real experimental conditions and filter func- 
tions will be considered elsewhere, while now the above 
heuristic arguments are to be extended to finite n. These 
lead to the formulae (8)-(10) again, but all the quantities 
become n-dependent and a non-zero jj, n i appears in (10): 



3 - In q 2 /V<l 2 + a l > 9 < 9n,M 



9n,M — 3/7„ 



1/2+ (l/4+(a n7n /3) 2 ) 



i/2- 



1/2 



(11) 
(12) 



Vn, q = JnQ [arsh (q/a n ) - arsh (l/a„)] + fi n>1 q . (13) 

These formulae contain three functions of one variable, 
j n , a n and jJ. n ,i, unknown yet. The function jj, n i is small 



2 



and tends to zero at n — ► oo. The function j n is constant 
up to a small varying term and given by 

j n =7 + 2/^,1 . (14) 

It is sufficient now to have formulae for i and Hn,2 
in order to determine all n n>q . For derivation of scale 
dependence of intermittency exponents, quite different 
ideas are required which are considered in [6]. Here are 
formulae for 1 and Hn,2 without derivation: 

A«n,i = - log 2 9n,i /n , g n ,i = 1 + 6 (l - 2~ hn ) , 

b = D 1 + l-H2-26 , <S<1; (15) 

9n,2 = 92 + (I- 92) (G 2 /2) n + 2 ( fln|1 - 1) ; 

</ 2 = l/(2-G 2 ) , G 2 = 2^. (16) 

The quantity a n is expressed explicitly in the terms of 
j„, Hn,i and /i„ )2 by 

a n = 2 L4 n (2 - A n ) (2A n - 1)] 1/2 / (A n - l) , 

A n = exp [{Hn,2 ~ 2/U n ,l) / (27n)] • (17) 

The singularity exponents a n , q = iJ,' n q &t n = I; ca and 
q = 0; g n ,M satisfy the relations 

2 3-«i, _ 2 «i,M + 7 ; 2 3_a ° = 2" M + 7 - 26 . (18) 

These formulae allow one to express all n n>q in the terms 
of the single parameter 7. 

• COMPARISON WITH EXPERIMENTS 
AND PHENOMENOLOGICAL MODELS: The 

experimental intermittency exponents n n>q for q ranged 
from 2 to 10 and the scale ratios x/y from 1.5 to 32 
(which corresponds to n ranged from 0.6 to 5) are pre- 
sented in paper [2]. The best fit to experimental Hn,2 
is reached at 7 = 0.296 . The corresponding values of 
the most important other parameters, calculated from 
the above formulae, are: ^2 = 0.214; 6 = 0.023; b = 
3.623; a M = 0.817 « 4/5; a = -0.124 « -1/8 . 

The theoretical intermittency exponents n n>q as func- 
tions of x/y = 2" at q = 1, 2, 10 are presented in the 
Fig.l-top (the order is clear as n n>q is increasing function 
of q). The circles are the experimental \i n q from Fig. 7 of 
[2]- 

The theoretical intermittency exponents n n>q as func- 
tions of q for x/y = 2" = 1, 1.5, 2, 4, 8, 32, 00 are pre- 
sented in the Fig.l-bottom by solid lines (the order is 
clear as n n>q is increasing function of n). The circles 
are the experimental n n>q from Fig. 7 of [2]. The pluses 
are \i q = q — (3^, where ( p are the experiment-based 
velocity structure functions exponents from the Fig. 3 of 
[3]. The 'x-marks' are experimental \i q from [4] (see be- 
low concerning the negative q). The dash-dot line is the 
Meneveau-Sreenivasan (MS) phenomenological model [4] 
corresponding in current notations to = 2, = 1/2 



(parameter is slightly adjusted, so that ^2 = 0.228, 
rather ^2 = 0.25 as in [4]). The dashed line is the She- 
Leveque (SL) phenomenological model [5] based on other 
set of experiments and corresponding in current notations 
to Doo = 1, = 2/3. 
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FIG. 1. Intermittency exponents fi n ,q as functions of 2" - 
top, and of q - bottom (see text for details). 

As is clear from the Fig.l-bottom, the limit n = 00 of 
the current theory is practically indistinguishable from 
both the SL and adjusted MS models in the range < 
q < 4. This range includes the whole range of data (ap- 
proximately < q < 3) that was used to justify SL- 
model. Noteworthy, that the data of [2] differ substan- 
tially from the x/y = oo-limit even at x/y = 32, where, 
for instance, Hn,2 ~ 0.16. There is no contradiction how- 
ever according to the current theory, as the finite size 
corrections \i q — pt nA decrease very slowly at n — ► 00. 

The theoretical dimensionalities D n>q as functions of 
q for x/y = 2 n = 1, 1.5, 2, 4, 8, 32, 00 are presented in 
the Fig. 2 by solid lines (the order is clear as D n>q is a 
decreasing function of n). The dash-dot and dashed lines 
represent the dimensionality D q defined as above and 
calculated according the formulae for \i q of the MS and 
SL models respectively. 
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FIG. 2. Dimensionalities D n<q as functions of g 

• NEGATIVE MOMENTS: The negative q could 
be considered in a similar manner as positive q, if one as- 
sumes that domains of low dissipation are associated with 
the same singular structures as hot spots. Such an as- 
sumption is justified for sufficiently anisotropic collapse, 
for which formation of hot spots includes substantial an- 
gular redistribution of local energy flux e near a singu- 
larity, and thus produces a "cold spot" right near a hot 
spot. Then the number of cold spots of a given intensity 
should be the same as the number of hot spots of the 
corresponding intensity. The selection of colder spots by 
taking higher negative moments q should be a uniform 
process in q, similar to that considered above for positive 
q. As a result similar formulae arise. If one assumes that 
the first two derivatives of D nq over q 2 are continuous at 
q = 0, it would imply that D n - q = D n q , i.e. D n q is an 
even function of q. 

Formula (6) for n n>q is applicable for q < as well. For 
D n ,q even in q, the first item in (6) is also an even function 
of q. Then the odd part of n n>q is a linear function of q, 
specifically, n n>q — n n ,-q = 2a n) o§ • Since a n o < 0, the 
intermittency exponents for a negative q are larger than 
for the positive q of the same absolute value. 

Thus extended to negative q the dimensionality and 
intermittency exponents are presented in Fig. 2 and Fig.l- 
bottom, as well as the extensions of MS and SL models. 
The SL model does not pretend to describe negative q 
and it is not related to experimental data in this region. 
The MS model is desired to fit experimental \i q for all q 
including negative ones. As is seen from Fig.l-bottom, 
the theoretical \i q agrees with the MS experimental data 
of [4] for negative q as good as the MS model. 

• FEW REMARKS: It is important to determine 
the experimental n n>q with a small step in q (which is 
a problem of the data processing, rather than the mea- 
surements). This would allow one to make an accurate 
differentiation and find experimental D n q . The accuracy 
requirements are clear from the Fig. 2 where the MS and 
my D q strongly diverge at negative q, although the MS 



and my \i q are close to each other in the same region of 
Fig.l-bottom. 

The first intermittency exponent \i n> \ is small, but it is 
important for the quantitative comparison between the 
theory and experiment. To the best of my knowledge, 
there is no direct experimental data on This is 

probably because the first intermittency exponent is usu- 
ally assumed to be zero. It is desirable to determine the 
experimental \i n> \ and to compare it with the above pre- 
dictions. 

There is a fundamental problem of applying filters 
more general than one-dimensional cuts currently used 
in nearly all experiments. 

When the intermittency exponents are known for all q, 
one can restore the probability distribution of the break- 
down coefficients M(r,x,y) = e(r, y)/e(r, x). It has a 
symmetry corresponding to the considered above hidden 
symmetry of the intermittency exponents. 

There is a connection between the scale-dependent in- 
termittency exponents and the velocity structure func- 
tions in developed turbulence. 
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